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^ Abstract 

Oh- Using gauge invariant quark Green's functions, defined with path-ordered gluon field 
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phase factors along polygonal lines, and functional relations among them, two compatible 

bound state equations of the Dirac type are established for quark-antiquark systems, 

each relative to the quark or to the antiquark of the system. The kernels of the bound 

state equations are defined through series of Wilson loop averages along closed polygonal 

contours and their functional derivatives on them. A sufficient criterion for spontaneous 

chiral symmetry breaking is derived, relating the Goldstone boson wave function in the 

zero total momentum limit with the scalar part of the gauge invariant quark two-point 



Green's function. 
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1 Introduction 

Gauge invariant quark Green's functions (GIQGF) [H [2], together with Wilson loops 
[21 IH El El El El E], represent the natural tools for the investigation of the properties 
of observable quantities in QCD [TQ1 [HI [121 EH lU EE51 dS|- Approaches using these 
ingredients meet, however, difficulties arising mainly from the nonlocal structure of the 
GIQGFs. 

In this respect, to optimize the methods of investigation based on GIQGFs, an ap- 
proach was undertaken by the present author with the aim of obtaining integrodifferen- 
tial equations that the latter would satisfy, in a parallel way as for the Dyson-Schwinger 
equations in the case of ordinary Green's functions [T71 [18] . This was possible in the case 
of two-point GIQGFs (2PGIQGF) defined with a path-ordered gluon field phase factor 
along polygonal lines between the quark fields [TS]. The 2PGIQGFs can then be clas- 
sified according to the number of segments their phase factor line contains. Functional 
relations are then obtained among these 2PGIQFs, which, together with the equations 
of motion relative to the quark fields, lead to an integrodifferential equation satisfied by 
the 2PGIQGF defined with one straight line segment for the phase factor. The kernel of 
this equation involves, with increasing complexity, a series of Wilson loop averages along 
polygonal contours and their functional derivatives. 

As a first step for the resolution of the above equation and the determination of 
the most important piece of the kernel, the case of two-dimensional QCD in the large- 
N c limit [20] was considered. The equation could then be solved exactly and analytically, 
displaying the main features of the spectral properties of the quark fields [21] : The quarks 
contribute to the 2PGIQGF like physical particles with positive energies, respecting the 
causality property; the singularities of the GIQGF are located on the positive real axis of 
the momentum squared variable (timelike region) and are represented by an infinite series 
of branch cuts. Although results obtained in two-dimensional theories cannot straight- 
forwardly be transposed into four dimensions, they underline here two features: i) The 
equation obtained for the 2PGIQGF is neither empty, nor unresolvable. With plausible 
assumptions about the properties of Wilson loops, it might also be analyzed in four- 
dimensions. In particular, the only part of the kernel that survives in two dimensions 
is precisely the simplest Wilson loop, corresponding to a triangular contour, ii) The 
resolution of the equation has provided new results, not known previously from more 
conventional approaches. 

The aim of the present paper is to enlarge the scope of investigations of the GIQGFs 
by also including in it the four-point GIQGFs (4PGIQGF), which allows us to study the 
bound state problem of quark-antiquark systems with a gauge invariant formalism. This is 
done by using again the functional relationships between GIQGFs with different numbers 



of segments on their polygonal lines. One then ends up with two bound state equations 
of the Dirac type, each relative to the quark or to the antiquark of the system. The two 
equations are compatible among themselves due to the validity of the Bianchi identities 
satisfied by the gluon fields. The kernels of these equations involve, as in the 2PGIQGF 
case, series of Wilson loop averages along polygonal contours and their derivatives, as well 
as the 2PGIQGFs of each quark field. A functional relationship is established between 
these kernels and those of the 2PGIQGFs. It is assumed, for later investigations, that the 
bound state wave functions satisfy usual spectral properties, based on the positivity of the 
energies of the quark and gluon fields and on causality, leading here to a generalization 
of the Deser-Gilbert-Sudarshan representation |22j . 

The bound state equations thus obtained allow us to further investigate the question 
of spontaneous chiral symmetry breaking. It is shown that in the chiral limit (massless 
quarks), the bound state equations possess a massless solution with zero total momentum 
if the 2PGIQGF possesses a normalizable nonvanishing scalar part in this limit. This 
result is the analogue of the one established by Baker, Johnson and Lee [23] for the 
Bethe-Salpeter equation [2U [25], [26] , and provides, prior to the resolution of the bound 
state equations, a sufficient criterion for chiral symmetry breaking in the case of bound 
states made of quarks with different flavors. 

The plan of the paper is the following. In Sec. [2] properties of the 2PGIQGFs are 
summarized. In Sec. [3] 4PGIQGFs are introduced and their properties are displayed. In 
Sec. HI the bound state equations are established. Sec. [5] deals with the question of the 
spectral representation of the wave functions. In Sec. El a criterion for chiral symmetry 
breaking is derived. Summary and concluding remarks follow in Sec. [71 An appendix is 
devoted to the normalization condition of the wave functions. 

2 Two-point Green's functions 

We summarize in this section the main results obtained for the 2PGIQGFs in Ref. |19j . 
We shall mainly be interested by path-ordered phase factors that are defined along polyg- 
onal lines in space (skew-polygonal), themselves being made of junctions of straight line 
segments. Designating by U(y, x) a phase factor along the straight line segment xy, with 
an orientation from x to y, a displacement of one end of the segment, while the other one 
is fixed, generates a displacement of all points of the segment with appropriate weight 
factors. We may characterize this as representing a rigid path displacement. Parametriz- 
ing linearly the segment with a parameter A, < A < 1, such that a point of the segment 
is represented as z(X), with z(0) = x and z(l) = y, a displacement of one end point of 
the segment gives rise to two types of contribution, the first coming from the end point 



itself and the second coming from the inner points of the segment. One has for the rigid 
path derivatives the formulas 

dU Q ya X) = -igA a {y)U{y, x) + ig(y - xf £ dX X U(y, z(X))F 0a (z(X))U(z(X),x), (2.1) 

dU ^ a X) = +igU{y, x)A a {x) + ig(y - xf £ dX (1 - A) U(y, z(\))F 0a (z(X))U(z(X),x), 

(2.2) 
where A is the gluon potential, F its field strength and g the coupling constant. In the 
above equations, the integrals represent the inner contributions of the segment. When 
dealing with gauge invariant quantities, the end point contributions are usually cancelled 
by similar contributions coming from neighboring segments or fields and it is the inner 
contributions of the segments that remain. We adopt for them the following notations: 

6 J^1 = ig{y _ x f £ dX X U(y, z(X))F Pa (z(X))U(z(X),x), (2.3) 

5 J^1 = ig{y _ x f £ dX (1 - A) U(y, z(X))Fp a (z(X))U(z(X),x). (2.4) 

Taking into account the orientation on U, the superscript + or — of the derivative variable 
indicates the segment on which it is acting when we are in the presence of two joined 
segments. Thus if we have the expression U(y,u)U(u,x), then the operator 5/5u + will 
act on U(u, x) only, through the end point u of the segment xu, while the operator S/Su~ 
will act on U(y,u). 

The vacuum expectation value of a Wilson loop (or, equivalently, the Wilson loop 
average) along a closed polygonal contour with n sides and n junction points X\, x 2 , ■ ■ ., 
x n , will be denoted W n and will be represented as the exponential of a functional F n 018]: 

W n = W(x n , x n _i, . . . , xi) = e Fn ( Xn > x »-i> " " " ' Xl ) = e Fn , (2.5) 

the orientation of the contour going from x± to x n through x 2 ; X3, etc. Then, the notation 
8F n /8x^ means that the derivative acts on the internal part of the segment XiXi + i with 
Xi + \ held fixed (x n+ i = xi), while SF n /Sxf means that the derivative acts on the internal 
part of the segment Xi_\Xi with %i-\ held fixed (x = x n ). 

The 2PGIQGFs with phase factors along polygonal lines can be classified according 
to the number of segments they contain. The 2PGIQGF with a phase factor line with n 
sides and n — 1 junction points ti, t 2 , . . ., £ n _i between the segments, is defined as 

S^) (2, x'; £ n _i,..., £1) = -— {ijj{x')U(x' ,t n -x)U{t n ^ 1 ,t n -- 1 ) ■ ■ ■U(t 1 ,x)ij){x)), (2.6) 

the ^s being the quark fields, with mass term m, belonging to the fundamental repre- 
sentation of the color gauge group SU(N C ) and the vacuum averaging being defined in 



the path integral formalism. (Spinor indices are not written and the color indices are 
implicitly summed.) The orientation of the path in S( n )(x, x'\ £ n _i, . . . , t\) runs from x to 
x', passing by t x , t 2 , . . ., t n _ x . 

The simplest 2PGIQGF corresponds to the case where n — 1, for which the points x 
and x' are joined by a single straight line: 
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S(i)(x,x') = S(x,x') = --^- (i/;(x')U(x',x)ip(x)). 



(2.7) 



(We shall generally omit the index 1 from that function.) A graphical representation of 
the 2PGIQGFs S^ and S^ is shown in Fig. [TJ 
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Figure 1: Graphical representation of the 2PGIQGFs 5(i) and S^. The solid lines rep- 
resent the quark field contractions, the dotted lines the phase factor along the polygonal 
lines and the arrows the orientation on them. 



For the internal parts of rigid path derivatives, we have the definitions 
8S( n ) (x, x'\ t„_i, ...,h) _ 



5x^~ 



-i- l^{x') U(x', t n _ 1 )f/(t„_ 1 , t n _ 2 ) . . . SU j h ;_ x) i>(x)), (2.1 



N, 



5x^- 



5SU^t^...M) = _j_ {W) 8U(£U-i) u{u _ u tn _ 2) u{tux) ^ x)y (2 0) 
5(i) and S( n ) satisfy the following equations of motion: 



(ij.d( x ) — m)S(x, x') = iS 4 (x — x') + iY 



SS(x,x') 
_ 5x»- ' 

S(x,x')(-ij. dw -m) = i5\x - x') - i 6 -S^p- 1 ^ 

(ij.d {x) - m)S {n) (x, x'- t n _i, . . . , t x ) = i5\x - x')e F ™( x ' *«-i> ■ ■ ■ ' **) 

+l1 Ix^l ' 

S (n) (x, x'- i n _i, . . . , tJi-ij. dp) -m) = i5\x - x')e F ^ x ' *«-i> • " • ' *0 



. 5S(n)(x,x'\t n - 1 ,. ■ ■ ,ti) 



7' 



(2.10) 
(2.11) 



(2.12) 



(2.13) 



The S7 n )S also satisfy equations related to the junction points t\, ..., t n _i on the 
polygonal line or more generally to local deformations of the paths. These involve the 
gluon field equations of motion and lead to equations related to the properties of phase 
factors and Wilson loops [U EJ El [7j . They should mainly be used for the determination of 
the expressions of the Wilson loop averages. In the present paper, the latter are assumed 
to be known and therefore the corresponding equations will not be considered. 

Multiplying the equations of motion (j2.12p and (12.131) with S(ti,x) and S(x',t n -i), 
respectively, and integrating with respect to x or x', one can establish functional relations 
between the various 2PGIQGFs. For S( n ), one has 

S (n) (x, x'; t n . u • • • , ti) = S(x, x') e^+i^'' *"-i. ■■■ i h,x) 

,6S(x,yi) , , 5 \ , , . 

+ { 7 51+ + S ( x i yi)j^^J 7 S (n+l)(yi, % I *n-l, ■■■,t 1 ,X) 

= S(x x') e^ n+1 ( x ' ^n-u ■ ■ ■ ' *1' x ) 
-S [n+1) (x, z i; x', t n . u ...,h) 7 ft ( ^h^ + T^IT s &> x ')) ■ ( 2 - 14 ) 

OZ-^ Z-^ 

[Integrations on intermediate variables are implicit and will not be written throughout 
this paper. Here, y\ and Z\ are integration variables.] By iterating these equations with 
respect to the higher S( n )8 of the right-hand sides and assuming that the last term rejected 
to infinity tends to zero, one ends up with a series expansion of any Sua (n > 1) in terms 
of S and its derivative and derivatives of logarithms of Wilson loop averages. This shows 
that among the infinite set of 2PGIQGFs with polygonal lines, only the first one, with 
one single straight line, is a genuine dynamical independent quantity; all others are in 
principle calculable from it, provided one knows to evaluate the rigid path derivative of 
S and the Wilson loop averages. 

The calculation of S proceeds from the equation of motion (j2.10p (or (12.111) ). It is 
then necessary to devise a method for evaluating the rigid path derivative 5S(x, x')/8x^~ . 
This is done by applying the rigid path derivative operator on both sides of Eq. (12. 14ft 
and repeating iteration operations. Specializing the result for S, one finds 



71=3 



/6S (y n - 3 ,y n - 2 ) 

to 2+ ' — 8yZ 



x '" x V — T^^+ — + S{y n -^yn-2) - a„_ 2 - j7 



OLn-2 



x f - f a „_i+ S{y n - 2 ,y n -i)l S( n) {y n - 1 ,x;x,y 1 ,...,y n -2 J 



(2.15) 



The right-hand side involves a series of terms in which the nth-order one contains 
S7 n ) and a Wilson loop with a polygonal contour with {n + 1) sides. One notices from 
the locations of the second-order derivatives acting on the Fs the absence of reducible 
type contributions in the corresponding expressions; the latter are expected to be part of 
the definition of the S7 n )S when expressed in terms of free propagators. The calculation 
should be completed by bringing all derivative operators to the right; the final form 
shows that the rath-order term contains globally n derivatives acting on the logarithm 
of the corresponding Wilson loop average and/or on the Green's function S (at most at 
first-order for the latter); furthermore, each derivative acting on the Wilson loop operates 
on a different segment from the others; this prevents the appearance of singularities 
arising from derivatives acting on the same point. The Wilson loop contributions have 
the characteristics of being irreducible and are classified into three categories: connected, 
crossed and nested [TU] . The form (I2.15P is the most convenient one for comparisons with 
other cases, as those of the 4PGIQGFs. 

Equation (12. 15[) can be considered as the analogue of the self-energy Dyson-Schwinger 
equation in the case of the ordinary Green's function [TTJ [18]. Defining the latter as 
S(x, x') = jf (ip(x) ij){x')), its equation of motion takes the form 

(ij-d(x) - m)S(x, x) = i5 4 (x - x') + E(x, y) S(y, x'), (2.16) 

where E defines the self-energy; it is a functional of the Green's function S itself, together 
with other Green's functions of interest, like that of the gluon field or of the photon field 
in the case of QED. The expression of T,[S] in terms of S and the other two-point Green's 
functions defines the Dyson-Schwinger equation. 

In the present case, however, one meets a more complicated situation in two respects. 
First, Eq. ( J2.15D involves in its right-hand side the whole set of 2PGIQGFs defined 
along polygonal lines, although all of them are ultimately expressible in terms of S. This 
feature, on the other hand, is an indication that the set of 2PGIQGFs along polygonal 
lines is closed, since no other types of contour are needed to reach the final equation for 
S. Second, the integrals that are present are not of the convolution type; they overlap 
all terms that accompany them; this is due to the presence of the Wilson loops, whose 
contours pass by all points that are present. Taking into account these facts, one needs to 
introduce matrix type self-energy operators, S mn , where the first index refers to the initial 
Green's function that is considered (S'(m)) and the second one to the Green's function S( n ) 
that appears in the right-hand side of the equation. With this definition, the equation of 
motion ( 12. lQj) can be schematically written in the form 

oo 

(i 7 .d {x) - m)S(x,x') = i5\x - a/) + £ {^n[S] * S {n) ), (2.17) 

ra=l 



where the star operation represents the integrals involved in the term containing Sr n ) in 
the right-hand side of Eq. (12.151) and the functional expression of EirjS'] is deduced from 
that equation by identification. 

The above results can also be applied with obvious transpositions to the equation of 
motion ( 12. lip . The evaluation of 5S(x, x')/5x M+ can be done in two ways. First, one 
might use the second expression of the functional relationship of Eq. ( I2.14p . In this 
case the iterative expansion is done in the reverse order to that of Eq. (I2.15p . The 
corresponding equation takes then the schematic form 



S(x,x')(-n. d&) -m) = i5\x - x') + £ (S {n) * Z ln [S}), (2.18) 

ra=l 

where EirJS] is deduced from Ei^S 1 ] by reversing the orders of appearance of various 
matrices and Green's functions and some of their arguments, changing the sign of terms 
containing an odd number of explicit 7-matrices and replacing the rigid path derivatives 
S/Sx^ and 5/5y" J by 5/5x'^ + and 8/8zj 3 , respectively. Second, one might use the 
same functional relations as for obtaining Eq. ( I2.15p . In this case, the only modification 
is the replacement of the the rigid path derivative <5/<5x M ~ by 5/5x M+ , while in Eq. (12. lip 
we have to take into account the global change of sign in front of SS(x, x')/5x M+ and the 
position of the matrix 7^ on the utmost right. We shall write the resulting expression of 
the equation in the form 



S(x,x')(-n. d&) -m) = i8\x -*') + £ (Xm[S] * S {n) ). (2.19) 

n=\ 

It is to be emphasized that once the operator EinfS 1 ] of Eq. (I2.17P (or its analogue of 
Eqs. ( I2.18p - ()2.19p ) has been evaluated by means of Eq. (I2.15p . then the Green's functions 
5(1), . . ., S(n), etc., have to be considered as ordinary complex functions of the variables x, 
x', yi, . . ., y n , etc., satisfying translation invariance and well defined Lorentz transforma- 
tion properties. The information contained in their phase factors along the rigid straight 
line segments and polygonal lines is now expressed by means of the corresponding Wilson 
loop averages and their rigid path derivatives, which, after evaluation, are themselves 
ordinary functions of their arguments x, x', etc. (the junction points of the segments). 
The straight line segments and the polygonal lines do no longer introduce additional de- 
grees of freedom, since their geometry is completely determined by the knowledge of the 
positions of the junction points of the segments. In particular, the operator d/dx of the 
Dirac operator in Eq. (I2.17P acts as an ordinary derivative operator on S. This is why, 
Eq. ( I2.17P and its analogues have the status of integrodifferential equations. 

In the rest of the paper we shall consider systems involving two quarks with different 
flavors and generally with different masses. To distinguish their individual Green's func- 
tions we shall introduce an additional index for their notation. Thus, they will be denoted 
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Si,(n) and S 2) ( n ), respectively corresponding to quark 1 and quark 2. The simplest Green's 
functions (with one straight line) S (= 5(i)) will be denoted Si and S 2 . 



3 Four-point Green's functions 

We now consider two different quark fields, labelled with indices 1 and 2, respectively, 
with mass terms rri\ and 7712- Four-point GIQGFs are constructed by including gluon field 
phase factors between the quark and the antiquark fields. Considering a polygonal line 
made of n segments and another made of one segment, we define the 4PGIQGF G< n ) as 

G( n ) a /3,l3'a'{Xl, X 2 \ X 2 , X x ] t n _i,t n _ 2l . . . , ti) 

= -jT- {i ! 2 p{.x 2 )U{x 2} t n _ 1 ) . ..£/'(ti,a;i)^i a (a;i)V'i a /(xi)t/'(xi,x / 2 )^'(4)). ( 3A ) 

where a, (3, /3', a' are the spinor indices of the quark fields. 
The simplest such Green's function is G^y. 

G(i)(x u x 2 ;x' 2 ,x[) = -— {^ 2 (x 2 )U(x 2 ,x 1 )i)i{xi)i) l {x' 1 )U{x l l ,x' 2 )i) 2 {x' 2 )). (3.2) 

A graphical representation of G(\) and G(s) is shown in Fig. [2j 





G(i 



(i) 



r (3) 



Figure 2: Graphical representation of the 4PGIQGFs G(±) and G(3). Same conventions as 
in Fig. [TJ 



In the following, for the study of the bound state problem, the points x[ and x 2 will 
be sent to —00 in time; they will then disappear by factorization from the bound state 
equations; this is why the classification of the 4PGIQGFs hinges here only on the line 
between the points x\ and x 2 . 



The 4PGIQGFs satisfy the following equations of motion: 
(ij.di - mi)G(n)(xi, x 2 ; x' 2 ,x[; t n -i, ■■■,h) 

= i5 A (xi - x' 1 )S2,(n+l)(X2,X2]t n -i, . . . ,ti,Xi) 

+iY l ^-r^G {n) (x 1 ,x 2 ;x' 2 ,x' 1 ;t n _ 1 ,...,t 1 ) (3.3) 

ox\ 



lii 



G(n)(xi, x 2 ; x' 2 , x'^tn-t, . . . ,ti)(-ij. 82 -m 2 ) 



iS (x 2 -x 2 )Si t ( n+ i)(xi,x' 1 ;x 2 ,t n - 1 
. I 
5x 



ij^+G( n) (x 1 ,x 2 ;x 2 ,x' 1 ;t n - 1 ,...,t 1 ) t ^j v . (3.4) 



t„-lX2 



(di = d/dxi and d 2 = d/dx 2 .) The 7-matrices that act on G from the left, act on its first 
spinor index (a), while those acting from the right act on its second spinor index (/3), as 
defined in Eq. (13. ip . 

Multiplying Eqs. ( 13.31) and (13.41) with S\{b\,X\) and S 2 (x 2 , t n -i), respectively, and 
integrating, one obtains functional relations between different 4PGIQGFs: 

G( n )(xi, x 2 ; x 2 ,x l ; £ n _i, . . . ,ti) — Si(xi, x l )S 2 ^ n+2 )(x 2 , x 2 \ £ n _i, . . . ,ti, Xi,x 1 ) 

+ { 7 gl+ — +5 , i(xi,j/i)=- 5I3 j7 1 G^+^iy^ x 2 ; x 2 , x^t^^ . . . ,t u Xi) 

— <~>l,(n+2)(^l) x t ; x 2 , X 2 ,t n _i, . . . , ti)S 2 (x 2 , x 2 ) 

h fSS 2 (z u x 2 ) S 



-G (n+ i)(xi, z x \ x 2 , x[] x 2 , t n _i, . . . , ti) 7 A ( — I gp 2 + ^-g^S^zi, x 2 )). 



(3.5) 



4 Bound state equations 



In order to obtain bound state equations, one considers in the 4PGIQGFs the limit of large 
timelike separations between the set of points (x±,tx, . . . , t n -i,x 2 ) and the set (x 2 , x[). In 
this limit, the Green's functions can be saturated by a complete set of hadronic states, 
among which are single mesons representing bound states of quarks and antiquarks. To 
simplify the analysis, one may either consider the large- N c limit, in which case only single 
poles survive [201 I2Z] , or simply neglect inelasticity effects to ensure the stability of the 
bound state. By appropriate projection operations [25] or limiting procedures to the pole 
position in the total momentum space it is possible to select one particular bound state 
among the whole set of intermediate states. 

We shall refer to the selected bound state with its total four-momentum P only, dis- 
carding other quantum numbers, which will not play any role in the following. (We assume 
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that the bound state is nondegenerate.) The wave functions are classified according to 
the number of straight line segments existing on the polygonal lines of the phase factors 
(n,j = l,2,...,oo): 

1 



®{j)0'a'(x 2 , #i! tj-i, ■■■ ,h) 



SNr 



^K 



They become in the simplest cases n — 1 and j '• = 1 

1 



$(i)(xi,x 2 ) 
$(i)(ar / 2 ,ar / 1 ) 



fNr 



< 0\T(<lP 2 p( X2 )U(x 2 , t„-l) • • • C/(*l, Xi)Vla(xi))|P >, 

(4.1) 
< P\T{^ lal {x' l )U{x' 1 ,t^ 1 ) . . . U{t x ,x 2 )^{x' 2 ))\U > . 

(4.2) 

(4.3) 
(4.4) 



Wc 



< 0|T(^ 2 (x 2 )C/(x 2 ,xi)^i(xi))|P >, 



<P|T(^ 1 K)f/(x , 1 ,x 2 )^ 2 (x 2 ))|0>. 



(The dependence of the wave functions on the total four-momentum P of the bound state 
is omitted from their arguments for notational simplification.) We note that the above 
wave functions for all ns describe the same bound state, but differ in their expressions 
due to their differences in their contents with respect to the phase factor lines. 

Taking large timelike separations between initial and final coordinates as described 
above, or going in total momentum space to the pole position, and selecting one bound 
state, the equations of motion (I3.3l) - fl3.4p of the 4PGIQGFs are transformed into wave 
equations (their inhomogeneous parts not contributing to the bound states poles): 



(ry.di -mi)$(„)(xi,x 2 ;t n _i, . . . ,h) = -H7 M ^-^$(n)(xi,x 2 ; t n _i, 

OX-^ 



^(n)(xi,x 2 ;t n _ x , . . . ,t!)(-ij. <9 2 -m 2 ) 



which become for the case n 



Sx 



®(n)(xi,x 2 ;t n . 



H^(n) 



1, . . . 



*1 



*l) 



Xltl 



t„-lX 2 
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(ij.di -m 1 )$ (1) (xi,x 2 ) = +z'7 /1 j-^ Z :$(i)(xi,x 2 ), 

OX^ 



$ ( i)(xi,x 2 )(-z7. <9 2 -m 2 ) 



'5x u o + 



$(i)(xi,x 2 )Y- 



(4.5) 

(4.6) 

(4.7) 
(4.8) 



A wave function $( n ) thus satisfies two independent Dirac type equations. They 
should, however, be compatible among themselves in order not to give rise to a vanishing 
solution. The compatibility condition is obtained by making the two Dirac operators act 
on the wave function in different orders and subtracting the corresponding results from 
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each other. Since the two Dirac operators commute among themselves, the result should 
be zero. One finds 

^W'WW^""^ - (4J) 

For n > 2, the commutativity of the two rigid path derivatives results from the fact 
that they operate on different segments in an uncorrelated way. For n — 1, since they 
operate on the same segment, they may act on the same point, giving rise to additional 
singularities. It can, however, be shown that because of the Bianchi identities satisfied 
by the gluon fields, even in this case the two operators commute [TO]. The two wave 
equations (I4.5|) - (l4.6p . or (I4.7|) - (l4.8p . are therefore compatible among themselves. 

In order to evaluate the interaction part of the wave equations, it is necessary to 
express the rigid path derivatives in terms of calculable kernels. We follow here a method 
similar to that used for the 2PGIQGFs (Sec. [2]). By selecting in Eqs. ( 13.51) the bound 
state contribution, one obtains the following two equivalent equations, expressing a wave 
function $( n ) in terms of $( n +i) and the 2PGIQGFs Si or S 2 : 

<&( n )(xi,x 2 ;t n - 1 ,. . . ,ti) = 

rn / \ c 

= -&( n +l)(xi,Zi]X2,t n -i,...,tx)'Y Pl ( 1 ^_ 2 + y^S 2 (z 1 ,X 2 )). (4.10) 

OZ^ OZ\ 

However, contrary to the 2PGIQGF case [Eqs. ( 12.141) ]. they do not contain the lowest- 
index wave function $m which could generate an iterative series and allow for the cal- 
culation of the $( n )S in terms of $(i)- The difficulty can be overcome by adding to the 
above equations identities involving $(i). Considering the equations of motion ( 12.121) and 
( I2.13p . multiplying them with $(i)(ti,a?i) and &(i)(x 2 ,t n -i), respectively, and integrating, 
one obtains the two equations 

$(i)(2M, x 2 ) e Fn +^ x ^ *"-i' •••'*!' Xl ) 

+ ( — j Ql + + $(1)^1,2/1)^) 7 Ql S 2 ,( n+ i)(yi,x 2 ;t n -i, ■ ■ .,h,xi) = 0, 

(4.11) 

$(1) (Xi, X 2 ) e Fn + 1 ^ *«-!' •■■>*!» X l) 



-Si,(„+i)(xi, 2 X ; x 2 , t n _i, . . . , ii) 7 ft ( — ( i } g^ + ^-^$(1) (21, x 2 )) = 0. 



j ^$(1)^1, x 2 ) 

(4.12) 
These can now be added, respectively, to the two expressions of $(„) in Eqs. ( 14.101) . 
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yielding the following new functional relations: 

$(„)(si, x 2 ; t„-i, • • • , h) = $ (1) (x l7 x 2 ) e^ 1 ^ 2 ' *»-i' • • • ' *i' x 

+ ( — J Ql + + ^W^IjS/i) j-51^) 7 Ql ^2,(n+i)(yi,a:2; * n -i, • • ■ ,<i,xi) 

+ ( ^ai+ yi) + gl(^.yi)^p) T" 1 $(n+l)(yi,S2; tn-1, • • • ,*l,Sl), (4-13) 

$(«)(si, z 2 ; t n _ 1; . . . , ti) = $(l)(Xl, x 2 ) e^n+i^a, *n-l, • • • , *i, xO 



<— 



-5i i(n+ i)(xi, zi; x 2 , t n _i, . . . , tO 7 ( ^§jz + j-^+ *(i) (*i, ^2)) 

0,2^ 0,2^ 

-$ (n+1) (x l5 Zl ; x 2: t n _ 1; . . . , tO 7 ft ( ^^'^ + J ^+S 2 (z 1 , x 2 )) . (4.14) 

They can be used to express, through iterative calculations, $( n ) in terms of $(i), S±, S 2 
and Wilson loop averages. They parallel relations (12.141) of the 2PGIQGFs, but with the 
additional complication that the iteration should be carried out simultaneously in $( n +i) 
and S( n+ i). 

The final step consists in using expressions (I4.13p - p.14p to bring the right-hand sides 
of the wave equations (14.50 - (j4.8p into a form where the wave functions appear as acted 
on by kernels made of Wilson loop averages and their derivatives, as well as 2PGIQGFs. 
The procedure follows similar lines as those adopted for the 2PGIQGFs pj)]. Considering 
in particular the wave equation (14.70 . one finds 

M> (1) (xi,x 2 ) SF 2 (x 2 ,x 1 ) 
S 2 F 3 (x 2 ,xi,yi) 



L ®(i)(xi,yi) l ai S 2 ,(2)(yi, x 2 ; xi) + Si(xi, yi) 7 Ql $(2)(l/i, x 2 ; xi) 
r/5$ (1) (xi,y0 , * / \ $ \ ai 5 2 F A (x 2 ,xi,y 1 ,y 2 ) 

-{( - 5y?1+ +*u(zi,vi)^)>y - 6xrh?+ 

xS 2 (yi, y 2 ) 7° 2 S U3) (y 2 , x 2 ; x x ,yi) 
fSSi(xi,yi) 6 \ a J 2 F 4 (x 2 ,x u y u y 2 ) 

x [ $(i) (Vi, V2) 7° 2 S 2y{3) (y 2 , x 2 ; x 1} y x ) + Si(j/i, y 2 ) T 2 $(3) (j/a, z 2 ; zi, j/0 ] } 

-LL 1 ai + +Si(xi,yi)j-^)i x ••• 

i= 4r=l fyl <52/l 7 
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X ( La,+ +$ W(^-l>2/r)l^r)7 r ( ? ,.a r+1 + + & for , j/r+l) 1~5^: 



*#*+ ww " -"°"&y?->' V ^ ;r+ — — S y%i 
ar+1 t 6 S 2 (y 3-3^3-2) q ( x _£ \ a ,_ 2 

5 2 F i+ i(x 2 ,xi,?/i,...,y i _i) 
x g^-g aj _ 1+ S 2 (y j -2 1 Vj-i) 7 '- 1 S 2 ,(j){yj-i, x 2 ; x h y h ..., y i _ 2 ) 

j=4 V %1 %1 ' 

/5S'i(?/ j _3,?/ i _2) 



X 



+ 5 1 ( % _3,^ 2 ) I -^-^)7 aj - 
5 2 F., + i(x2,xi,2/i,...,2/j_i 






foi 5t/" J 



$(i)(%-2, 2/j-i) 7 aj_1 ^2,0)(%-i, x 2 ; xi, 2/1, . . . , %_ 2 ) 
+5 , i(%--2, %-i) 7 aj_1 $(,-) (j/j-i, ^2; Xi, j/i, . . . , j/j_ 2 ) 1 . (4.15) 



In the double sum with respect to r and j , yo means Xi ; when r = 1 , factors with S\ and 
derivatives do not exist on the left of $(iy, when r = j — 2, factors with £2 and derivatives 
do not exist on the right of &(i), except the one preceding S 2 uy 

The kernels appearing in the right-hand side of Eq. (I4.15P can be expressed through 
functional relationships in terms of those of the 2PGIQGFs. Going back to Eq. (12.17ft 
and to the definition of the operator Xi„[.S] there, to which we assign now a new index, 
1 or 2, according to its content of quark 1 or quark 2, the wave equation (14. 7p can be 
rewritten in the form 

(irdt - mi)$ (1) (x 1 , x 2 ) = £ (^, ln [Si] * *(„)) + E ( ^ Sl g N * * $ d)) * S Un)) , 

(4.16) 
where the double-star operator means that after functionally taking the derivative of 
Si,in[Si] with respect to Si, the latter is replaced in the same place by $(1) and all Sis 
on the right of $(1) are replaced by S 2 . 

A similar expression can also be derived for the right-hand side of Eq. (14. 8p . using 
either of the forms ( 12~Tgj) or (12TT3]) . 

Equation (I4.15P parallels Eq. (12. 15[) of the 2PGIQGF. High-index wave functions 
and 2PGIQGFs should in principle be eliminated in terms of the lowest-index ones using 
the functional relations (I4.13p - (l4.14p and (I2.14p . Derivative terms of $n) that appear 
in the right-hand side of Eq. (14.15ft in high-order terms could be eliminated using the 
first terms of the expansions through an iterative procedure. Expansion ( 14.15ft should 
be completed by bringing all derivative terms to the right; the result is very similar to 
what is obtained in the 2PGIQGF case. Connected parts of the Wilson loop averages of 
high-order terms involve an increasing number of derivatives, each on a different segment 
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of the corresponding polygonal contours 

In perturbation theory, each derivative of a Wilson loop introduces one multiplicative 
power of the coupling constant. Therefore, for the perturbative regime, terms with the 
smallest number of derivatives would be the dominant ones. For large-distances, Wilson 
loop averages are expected to be dominated by minimal surfaces [3J EJ [TU]. Here also, 
the dominant terms are those with the smallest number of derivatives. This suggests 
that the expansion (I4.15P has, at least formally, a perturbative structure, the decrease in 
magnitude of the terms being estimated by the global number of the derivative operators. 
The first term of the above expansions involves a Wilson loop with one derivative, which, 
however, vanishes for symmetry reasons. The first leading term of the expansion is then 
the two-derivative term, involving a Wilson loop along a triangular contour. In two- 
dimensional QCD in the large- N c limit, this term is actually the only one that survives 
in the above expansions and brings an indirect confirmation to the previous analysis. 

Our remark (in the paragraph following Eq. (12.191) ) concerning the interpretation of 
Eq. (12.171) and its analogues also applies to Eq. (I4.16P and its analogues. Once the rigid 
path derivative of the wave functions has been evaluated by means of Eq. (I4.15P or similar 
ones, then the wave functions have to be considered as ordinary complex functions of their 
arguments xi, £2, h, ■ ■ ■, t n , P, etc. The wave equations of the type of (14. 16ft become 
integrodifferential equations. 

For completeness, the normalization condition of the wave functions is presented in 
the appendix. 

5 Spectral representation 

A two-particle wave function of colorless fields is a vertex function and satisfies, on the 
basis of the properties of causality and positivity of physical particle energies, an inte- 
gral representation known as the Deser-Gilbert-Sudarshan (DGS) representation [22]. 
One might try to apply the DGS analysis to the presently defined wave functions (14. ip 
and (14.31) . However, an immediate difficulty arises from the fact that intermediate states 
needed for the determination of their singularities are necessarily colored objects. On the 
other hand, the equations satisfied by the GIQGFs do display momentum-space singular- 
ities generated by the free quark propagators present in them. The same difficulty also 
appeared in the 2PGIQGF case. 

To explain the presence of singularities in the Green's functions it is necessary to admit 
that completeness sums may be considered with colored quark and gluon states with pos- 
itive energies, irrespective of the fact that the latter may not be observable as asymptotic 
free states. It is the solutions of the corresponding equations which should ultimately de- 
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termine their precise properties. This hypothesis was applied to the 2PGIQGF case and 
led to a generalized Kallen-Lehmann representation [28j I2H1 [19] . The analytic resolution 
of the 2PGIQGF equation in two-dimensional QCD in the large- N c limit has confirmed 
that hypothesis: The quarks contribute to the spectral singularities like physical parti- 
cles, with positive energies, with the difference that their singularities are no longer simple 
poles, but an infinite number of branch points with a stronger power than poles [21] . 

It is therefore reasonable to continue to apply the same approach to the study of 
the spectral representation of the wave functions. Considering the wave function $m 
[Eq. ( 14. 3j) ]. one can repeat the same analysis as in Ref. pjj] for the 2PGIQGF. The 
presence of the gluon field phase factor introduces an infinite series of additional powers 
of the denominator of the dispersive integral. Defining total and relative coordinates as 
X = \{x\ + X2) and x = X\ — Xi, one can factorize the plane wave part of the wave 
function and introduce the internal wave function as 

$ w (P, Xl ,x 2 ) = e- %RX <t> { i){P,x). (5.1) 

0(i) (P,x) satisfies a generalized form of the DGS representation, which we write here, 
for simplicity, for the total spin-0 case, ignoring the spinorial content: 

The lower bound of the spectral variable s' may actually depend on the quark masses, on 
P 2 and on /3. The above general form might still be simplified through integrations by 
parts and recombinations. In particular, the sum might lead to a global fractional power. 
In two-dimensional QCD, the resulting power of the denominator for the 2PGIQGF case 
was found 3/2 |21j . 

Representation ( 15. 21) . or a simpler version of it, could be used for the resolution of the 
wave equations. 

6 Chiral symmetry breaking 

In the limit of Nf massless quarks, it is expected that the chiral SU(Nj) x SU(Nf) symme- 
try group of QCD undergoes a spontaneous breakdown to its diagonal subgroup SU(Nf) v 
[30] |3T| |32| [33] . On the other hand, the Nambu-Jona-Lasinio model [31] underlines the 
intimate relationship that exists between chiral symmetry breaking and the dynamical 
mass generation of fermions. Baker, Johnson and Lee considered the case of QED and 
showed that a similar relationship also exists there: the generation of a nonzero mass term 
in the fermion self-energy in the massless limit entails the existence of a zero-mass pseu- 
doscalar bound state solution of the Bethe-Salpeter equation [22] • They argued, however, 
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that such a solution might not necessarily correspond to an observable boson. The latter 
phenomenon is actually a consequence of the axial anomaly problem in abelian sectors 
[35| [36] and should be avoided by considering nonsinglet sectors of a nonabelian chiral 
group. 

In this section we aim at showing that the relationship between dynamical mass gen- 
eration of fermions and chiral symmetry breaking also exists in QCD. To this end, we 
consider the case of massless quarks 1 and 2; in this limit the labelling of the 2PGIQGFs 
with the quark indices 1 and 2 becomes irrelevant and we may discard them from our 
notations. The 2PGIQGF S^ is decomposed into vector and scalar parts [19]. The part 
that is most sensitive to chiral symmetry breaking is the scalar one. We isolate it by 
taking the anticommutator of S^ with the 75 matrix. It satisfies the equations 

ij.di[l5,S(i)(Xi,X2)] + = 27 M -— 3[7 5 , 1 S( 1) (Xi,X 2 )] + , (6.1) 

[7 5 ,5 (1) (xi,x 2 )] + (-z7. <9 2 ) = -«^Ti+[75,5 ( i ) (x 1 ,a; 2 )] + 7 /i , (6.2) 

where [, ] , represents the anticommutator. 

The right-hand sides of Eqs. (I6.ip - ()6.2p are obtained from Eq. (12.151) and its adjoint 
with respect to x 2 and more schematically from Eqs. fl 2 . 1 7D - ( 12 . 1 9 D . We have for the rigid 
path derivative with x\ 

I l o , m SF 2 (x 2l x 1 ) 

|^[75, 5 ( i)(xi, x 2 )] + = — jT-jz — [75, S {1) {xi, x 2 )\ + 

5 2 F 3 (x 2 ,x 1 ,y 1 ) r . 
5x^5 ai+ \^ S ^ Xl ' yi ^+ 7 * S 2,(2)(yi,a; 2 ;xi) 

+Si(x 1 ,yi)>y ai [75,5 , ( 2 )(2/i,x 2 ;xi)] + } +... , (6.3) 

where the dots correspond to similar terms coming from the contribution of the sum 
contained in Eq. (j2.15p . In schematic form, the equation relative to X\ is 

•Tftte, S(i)U = E (Sm[5 {1) ] * [75, %)]+) + f f( ^"^ (1)] * *[ 75 , S {1) } + ) * S (n) ) , 

n=l n=2 \ 0D (1) / 

(6.4) 
where Ei n and the star operation were introduced in Eq. (12.171) . while the definition of 
the double-star operator is the same as in Eq. (I4.16p . with the only difference that the 
quark indices, 1 or 2, are now removed from the Green's functions. 

Comparing the above equation with Eq. ( I4.16P (in the limits ni\ = m 2 = and 
removing the quark indices 1 or 2), one immediately concludes that [jb,S(x)(xi,x 2 )] + 
satisfies the same equation as $(i)(xi, x 2 ), provided the following correspondences are also 
done: [75, 5 , ( n )]+ - > ^(n)> n = 1, 2, . . . , 00, the general case of n being itself established 
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from the equations of Sr n y A similar conclusion also holds for the adjoint equations 
corresponding to Xi- On the other hand, the 2PGIQGFs do not depend on the bound 
state total momentum vector P. Therefore, the above correspondence is possible only if 
the $s are independent of P; this is possible if P = in the corresponding bound state, 
which in turn implies P 2 = 0. 

We thus arrive at the conclusion that if the 2PGIQGF 5(i) has, in the massless quark 
limit, a nontrivial normalizable scalar part, then the latter, multiplied by the 75 matrix, is 
a solution of the bound state equation with zero mass, in the limit of zero total momentum. 
This is an indication of the existence of a pseudoscalar Goldstone boson in the bound 
state spectrum. The complete expression of the corresponding wave function for nonzero 
P should be searched for from the bound state equations themselves; that part is not 
given by the 2PGIQGFs. 

It is worth noting that the scalar part of the 2PGIQGF Sm does not necessarily 
correspond to a conventional mass term which would give rise to a pole structure in the 
quark Green's function. Rather, the property that quarks are confined suggests that it 
would possess a more complicated structure. In two-dimensional QCD, the resolution 
of the problem led to the appearance of an infinite set of branch points, at dynamically 
generated mass values Mi, M 2 , . . ., M n , . . ., with stronger singularities than simple poles 

The above criterion for chiral symmetry breaking is only a sufficient one. This is 
a consequence of the fact that the equations satisfied by $ are linear in $ (given the 
expressions of the 5s), while the equations satisfied by [75, 5] + are nonlinear, since the 
latter are also contained in the 5s present in the equations. 

7 Summary and concluding remarks 

Phase factors along polygonal lines allow a simple classification of gauge invariant quark 
Green's functions according to the number of segments they contain on the lines and per- 
mit a systematic investigation of their properties through their equations of motion. The 
latter can be reexpressed as integrodifferential equations in which the kernels are essen- 
tially represented by Wilson loop averages along polygonal contours with an increasing 
number of sides and derivatives. 

This approach was applied to the cases of two-point and four-point Green's functions 
in QCD, leading in the latter case to bound state equations for quark-antiquark systems. 
The functional relationships between the kernels of the bound state equations and those 
of the two-point functions were displayed. A sufficient criterion for spontaneous chiral 
symmetry breaking was derived, relating the Goldstone boson wave function in the zero 



total momentum limit with the scalar part of the two-point Green's function. 

The idea of relating the perturbative degree of a kernel with the number of sides 
and derivatives of the Wilson loop contours offers promising perspectives for practical 
applications of the equations obtained thus far. In this case, the dominant terms of 
the interaction kernels would come from the simplest contours and the least number of 
derivatives. That feature is also manifest in two-dimensional QCD in the large- N c limit. 
An analysis of the structure of the present equations in two dimensions would provide us 
with a simplified framework for the understanding of various mechanisms at work and this 
in turn might serve as a helpful guidance for future resolutions of the relevant problems 
in four dimensions. 
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A Normalization condition 

The normalization condition of the wave functions is usually obtained by acting on the 
four-point Green's function with its inverse. In the present case, the inhomogeneous parts 
of the equations of motion (13.31) and ( 13. 4ft being two-point Green's functions, one needs 
to use, for Gm, simultaneously both equations. One obtains: 

i- 
A A 

(ry.di -m l - i'j'" T -—) G { i ) (x 1 ,x 2 ;x' 2 ,x[) (-17. 82 -m 2 -vf 



<fo£-' wv " " z ' i/v ' ' 5x v 2 + ' 

= i 2 5\ Xl - x\)5\x 2 - x' 2 ). (A.l) 

[F 2 (x 2 ,x 1 ) = l.] 

After calculating the effects of the rigid path derivatives, one passes to total momentum 
space P by taking the Fourier transform with respect to the total coordinate difference 
X — X' = \{x\ + x 2 ) — \{x' x + x' 2 ), while remaining in the relative coordinate spaces 
x = (x\ — x 2 ) and x' = (x[ — x 2 ). Designating by L(P,x) the operator acting in the 
left-hand side of Eq. (1A.1I) . the latter can be written in condensed form as 

L(P, x) * G {1) (P, x, x 1 ) = i 2 5 4 (x - x'), (A.2) 

where the star notation represents here the complete series of terms obtained by the 
action of L(P,x) on Gtu as in Eqs. (12.151) and (14.151) and which involves linearly the set 
of Green's functions Gt n ) (n = 1, 2, . . .). 
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We specify with the label k one of the bound states of the system and isolate its 
contribution and that of its antiparticle in Gm from the rest: 

G(i) (P, x, x ) = P 2_ P 2j rie + lR Wk (P^^), ( A.3) 

where the 0s represent the internal part of the wave functions, as defined in Eq. (15. ip . 
[The chronological product being evaluated with the time components of the xs, we have 
also Pk = P-] Decomposition ( 1A.3J) is replaced in Eq. ( 1A.2J1 and then both sides of the 
latter are multiplied with (f>m k (x) and integrated with respect to x, with the trace on the 
spinor indices taken. Taking into account the fact that the operator L(P k ) annihilates 
0(i)fc, one obtains 

— /L(P x) — L(P k x)\ — — 

<P(i)k(x)( p 2 _p2 + i g ) * ^(i)k(x) <P(i) k {x') + <f> {1)k {x)L(P,x) *R( 1)k (P,x,x') 

= ^(i)k(x). (A.4) 

At this stage, when working with conventional Green's functions, one takes the limit 
P 2 — > P k . The convolutive nature of L(P) then allows it also to act on the left and the 
remainder contribution disappears from the equation. In the present case, the operator 
L(P) acts specifically on the right and without further information it is not entitled to 
be converted to the left. 

To go further, we consider, as in Eq. ( lA.ip . equations of motion of Gm concerning its 
two right arguments and define correspondingly an operator L (P) acting from the right: 

G {1) (P, x", x)* L (P, x) = i 2 5\x" - x). (A.5) 

From Eqs. (IA.2J) and (1A.5I) one deduces 



G {1) (P,x",x)L(P,x) *G m (P,x,x') = G {1) (P,x",x)* L {P,x)G {x) {P,x,x') 

= i 2 G {l) (P,x",x'), (A.6) 

which entails a weak form of conversion of the operator L(P) from righ to left. Since this 
result is true for any P, x" and x', one might adopt the assumption that it remains true 
also for parts of Gn\. Adopting the latter assumption, and taking in Eq. (1A.4|) the limit 
P 2 — >■ P k , one obtains the two relations 



1 r 4 - dL{Px) 

-Jdx<P {t)k {x) gp2 



p2=Pi *<f>(i)k(x) = 1 i (A.7) 



J d 4 X(p {1)k (x)L(P k ,x) * R {l)k (P k ,x,x') =0, (A. 



which display the normalization and orthogonality conditions. In Eq. (1A.7I) . the derivative 
of L(P, x) with respect to P 2 is understood in the sense that one first evaluates the effect 
of L(P) on <f)(i)k and then takes in the resulting kernels the corresponding derivative. 



20 



Isolating in Rmk the contribution of another bound state, specified by a label m, 
and assuming its independence from the rest, one obtains the more precise orthogonality 
relation 

d 4 £0(i)fcO)( fc ' 2 _ p2 m ' — ) * 0(i)mO) = 0, m^k. (A.9) 

k m 



We notice that the free part of the normalization kernel in Eq. (1A.7J) . coming from 



the two Dirac operators, is the same as for the Bethe-Salpeter wave functions [26]. In 
the nonrelativistic limit, decomposing into 2x2 components, those which are dominant 
satisfy the properties 7 O = — </>7o = <$>■ 
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